Fourier Analysis Of Time Series Research Articles

An inflection point-based method for estimating metrics of mangrove phenology combining climatic factors and Landsat NDVI time series

Abstract The present research evaluated the prospects of utilizing rainfall and temperature combined with Landsat-8 derived HANTS (Harmonic ANalysis of Time Series) reconstructed NDVI for estimating the metrics of the mangrove phenology. The selected period of the study was from 2013 to 2020 for the Pichavaram mangroves of Tamil Nadu. The NDVI and ERA5 (ECMWF Re-Analysis) datasets of rainfall and temperature were the input datasets for developing the new algorithm. The ‘z-score sum’ provided a measure of the cumulative impact of rainfall and temperature, displaying its most negative value coinciding with the peak positive value of the NDVI time series datasets. The algorithm developed for phenological metrics estimation identified the common inflection points of the z-score sum and NDVI curves. The temporal analysis of metrics revealed the average Length of Season (LoS) as 230 days. The metrics also identified the drought year 2016 with the shortest LoS and the least Gross Primary Productivity (GPP) values. The analysis showed the influences of the preceding year’s monsoon rainfall on the GPP values of the later part of the phenological cycle. The temperatures during the days of PoS were found to be the optimum temperature for the growth of mangroves.

Estimating the CSLE Biological Conservation Measures’ B-Factor Using Google Earth’s Engine

The biological conservation measures factor (B) in the Chinese Soil loss Equation (CSLE) model is one of the main components in evaluating soil erosion, and the accurate calculation of the B-factor at the regional scale is fundamental in predicting regional soil erosion and dynamic changes. In this study, we developed an optimal computational procedure for estimating and mapping the B-factor in the Google Earth Engine (GEE) cloud computing environment using multiple data sources through data suitability assessment and image fusion. Taking the Yanhe River Basin in the Loess Plateau of China as an example, we evaluated the availability of daily precipitation data (CHIRPS, ERA5, and PERSIANN-CDR data) against the data at national meteorological stations. We estimated the B-factor from Sentinel-2 data and proposed a new method, namely the trend migration method, to patch the missing values in Sentinel-2 data using three other remote sensing data (MOD09GA, Landsat 7, and Landsat 8). We then calculated and mapped the B-factor in the Yanhe River Basin based on rainfall erosivity, vegetation coverage, and land use types. The results show that the ERA5 precipitation dataset outperforms the CHIRPS and PERSIANN-CDR data in estimating rainfall and rainfall erosivity, and it can be utilized as an alternative data source for meteorological stations in soil erosion modeling. Compared to the harmonic analysis of time series (HANTS), the trend migration method proposed in this study is more suitable for patching the missing parts of Sentinel-2 data. The restored high-resolution Sentinel-2 data fit nicely with the 10 m resolution land use data, enhancing the B-factor calculation accuracy at local and region scales. The B-factor computation procedure developed in this study is applicable to various river basin and regional scales for soil erosion monitoring.

Investigating the capability of the Harmonic Analysis of Time Series (HANTS) algorithm in reconstructing time series images of daytime and nighttime land surface temperature from the MODIS sensor

Investigating the capability of the Harmonic Analysis of Time Series (HANTS) algorithm in reconstructing time series images of daytime and nighttime land surface temperature from the MODIS sensor

Proof of evidence of changes in global terrestrial biomes using historic and recent NDVI time series

Climate change affects plant dynamics and functioning of terrestrial ecosystems. This study aims to investigate temporal changes in global vegetation coverage and biomes during the past three decades. We compared historic annual NDVI time series (1982, 1983, 1984 and 1985) with recent ones (2015, 2016, 2017 and 2018), captured from NOAA-AVHRR satellite observations. To correct the NDVI time series for missing data and outliers, we applied the Harmonic Analysis of Time Series (HANTS) algorithm. The NDVI time series were decomposed in their significant amplitude and phase given their periodic fluctuation, except for ever green vegetation. Our findings show that the average NDVI values in most biomes have increased significantly (F-value<0.01) by 0.05 ndvi units over during the past three decades, except in tundra, and deserts and xeric shrublands. The highest rates of change in the harmonic components were observed in the northern hemisphere, mainly above 30° latitude. Worldwide, the mean annual phase reduced by 9° corresponding to a 9 days shift in the beginning of the growing season. Annual phases in the recent time series reduced significantly as compared to the historic time series in the five major global biomes: by 14.1, 14.8, 10.6, 9.5, and 22.8 days in boreal forests/taiga; Mediterranean forests, woodlands, and scrubs; temperate conifer forests; temperate grasslands, savannas, and shrublands; and deserts, and xeric shrublands, respectively. In tropical and subtropical biomes, however, changes in the annual phase of vegetation coverage were not statistically significant. The decrease in the level of phases and acceleration of growth and changes in plant phenology indicate the increase in temperature and climate changes of the planet.

Long-term monitoring and phenological analysis of submerged aquatic vegetation in a shallow lake using time-series imagery

Submerged aquatic vegetation (SAV) is a key functional group for the restoration of wetlands and aquatic ecosystems. It plays a vital role in the stability of ecosystem structure and function in shallow lakes and reservoirs. With Baiyangdian Lake as the research area, Landsat TM/ETM+/OLI and Sentinel-2 imagery datasets were used to explore the phenological characteristics of SAV over different time scales and construct a phenological map of the SAV groups. The long-term spatial distribution characteristics and trends of SAV from 1986 to 2021 were analyzed. The results show the following: (1) the harmonic analysis of time series (HANTS) eliminated abnormal observations and noise. Normalized difference vegetation index (NDVI) time series curves generated after HANTS accurately reflected the phenological characteristics of SAV. The early spring submerged aquatic vegetation (SAV1), such as Potamogeton crispus, entered the germination period in late February at the earliest (DOY = 53). The mature period ran from April to May, after which SAV1 began to decline gradually. The length of the growing season (LOS) was about 110 days. SAV2, which was represented by Ceratophyllum demersum, Myriophyllum verticillatum, and Potamogeton pectinatus, germinated in early April (DOY = 99). The mature period ran from July to September. The length of the growing season was 210 days. (2) The decision tree constructed based on the NDVI and modified normalized difference water index (MNDWI) identified SAV with an accuracy of 89.7%. (3) The distribution range of the SAV changed dramatically during 1986–2021. According to the area change, it can be divided into a shrinkage period, expansion period, degeneration period, and initial period of recovery. The succession trend of SAV provide a theoretical scientific basis and technical support for the ecological restoration and management of shallow lake ecosystems.

Monitoring vegetation dynamics (2010–2020) in Shengnongjia Forestry District with cloud-removed MODIS NDVI series by a spatio-temporal reconstruction method

Shengnongjia Forestry District is the national natural reserve with abundant biological resources in China. Monitoring its variations of vegetation for the intimate connection with eco-environmental changes is of great significance. In this paper, the 16-day composite MODIS normalized difference vegetation index (NDVI) products (MOD13A1) with 500 m resolution from 2011 to 2020 were selected to investigate the vegetation dynamics in Shengnongjia Forestry District. To alleviate the cloud contamination in NDVI products, a spatio-temporal prefill method with harmonic analysis of time series (ST-HANTS) is proposed. ST-HANTS first prefills the raw NDVI time series using spatio-temporal information to reduce data gaps, which effectively improves the reconstruction performance of the subsequent HANTS algorithm. In the simulation experiments, ST-HANTS shows the highest average correlation coefficient and the lowest average root mean square error compared with other mainstream methods, including HANTS, Savitzky–Golay filter, wavelet transform, and data assimilation. The reconstruction curves are close to the upper envelope of the NDVI time series, which is more consistent with the vegetation phenology and can effectively capture the key points in the growth cycle. By analyzing the cloud-free NDVI time series reconstructed with ST-HANTS, results reveal the overall trend of Shennongjia vegetation coverage is high in the middle while low at the edge. Except for the population centers and Hongping Airport, the NDVI of most areas is above 0.7 and shows a remarkable increasing tendency. Moreover, the fluctuation degree of NDVI in the whole study area is low, indicating that the ecological environment of Shennongjia is relatively stable. The vegetation variations are influenced by land surface temperature and precipitation, and the vegetation growth response to precipitation exhibits an apparent hysteresis.

Gas hold-up analysis in an unsteady stirred vessel by means of infinite series

Abstract The use of a liquid level sensor made it possible to measure changes in gas hold-up over time in a stirred tank during unsteady mixing. These results were subjected to Fourier time series analysis and a model of gas hold-up changes in time was proposed. It allowed one to determine the model value of gas hold-up, which can be useful for characterizing gas hold-up during unsteady mixing and as a comparison to gas hold-up during steady mixing. The characteristic frequency was also determined, which corresponds to about twice the oscillation frequency. Model gas hold-up values for coalescing and non-coalescing systems were compared. Moreover, the change of the gas hold-up at constant maximum stirrer rotation frequency and variable gas flow rate for different oscillation frequencies was investigated.

A scalable software package for time series reconstruction of remote sensing datasets on the Google Earth Engine platform

ABSTRACT Spatiotemporal residual noise in terrestrial earth observation products, often caused by unfavorable atmospheric conditions, impedes their broad applications. Most users prefer to use gap-filled remote sensing products with time series reconstruction (TSR) algorithms. Applying currently available implementations of TSR to large-volume datasets is time-consuming and challenging for non-professional users with limited computation or storage resources. This study introduces a new open-source software package entitled ‘HANTS-GEE’ that implements a well-known and robust TSR algorithm, i.e. Harmonic ANalysis of Time Series (HANTS), on the Google Earth Engine (GEE) platform for scalable reconstruction of terrestrial earth observation data. Reconstruction tasks can be conducted on user-defined spatiotemporal extents when raw datasets are available on GEE. According to site-based and regional-based case evaluation, the new tool can effectively eliminate cloud contamination in the time series of earth observation data. Compared with traditional PC-based HANTS implementation, the HANTS-GEE provides quite consistent reconstruction results for most terrestrial vegetated sites. The HANTS-GEE can provide scalable reconstruction services with accelerated processing speed and reduced internet data transmission volume, promoting algorithm usage by much broader user communities. To our knowledge, the software package is the first tool to support full-stack TSR processing for popular open-access satellite sensors on cloud platforms.

Reconstructing NDVI and land surface temperature for cloud cover pixels of Landsat-8 images for assessing vegetation health index in the Northeast region of Thailand.

Critical applications of satellite data products include monitoring vegetation dynamics and assessing vegetation health conditions. Some indicators like normalized difference vegetation index (NDVI) and land surface temperature (LST) are used to assess the status of vegetation growth and health. But one of the major problems with passive remote sensing satellite data products is cloud and shadow cover that leads to data gaps in the images. The present study proposes temporal aggregation of images over a short time span and developing short span harmonic analysis of time series (SS-HANTS) and pixel-wise multiple linear regression (PMLR) algorithms for retrieving cloud contaminated NDVI and LST information from Landsat-8 (L8) data products, respectively. The developed algorithms were applied in the northeastern part of Thailand to recover the missing NDVI and LST values from time series L8 images acquired in 2018. The predicted NDVI and LST values at artificially clouded locations were compared with the corresponding clear pixel values. Additionally, the model predicted LST and NDVI values were also compared with MODIS LST and NDVI datasets. The calculated root mean square (RMSE) values were ranging from 0.03 to 0.11 and 1.50 to 2.98°C for NDVI and LST variables, respectively. The validation statistics show that these models can be satisfactorily applied to retrieve NDVI and LST values from cloud-contaminated pixels of L8 images. Furthermore, a vegetation health index (VHI) computed from cloud retrieved continuous NDVI and LST images at province level shows that most of the western provinces have healthy vegetation condition than other provinces in the northeast of Thailand.

Combing remote sensing information entropy and machine learning for ecological environment assessment of Hefei-Nanjing-Hangzhou region, China

Urban ecological environment is the basis of citizens’ survival and development. A rapid and objective urban ecological environment assessment (UEEA) plays an important role in the urban sustainable development and environment protection. This study established an improved urban ecological comfort index (UECIIMP), which is based on our previous UECI and fully composed of four remote sensing indicators: normalized difference vegetation index (NDVI), normalized difference built-up index (NDBI), land surface temperature (LST), and aerosol optical depth (AOD), representing the greenness, dryness, heat, and atmospheric turbidity, respectively. Combining the entropy method and random forest (RF) algorithm, the weights of four indicators were calculated. To improve the accuracy of UECIIMP, the gap-filled quarterly mean results of each indicator with 30m resolution were obtained using the harmonic analysis of time series (HANTS) method and spatial-temporal information fusion based on non-local means filter (STNLFFM). UECIIMP was applied to the Hefei-Nanjing-Hangzhou Region to explore its spatiotemporal changes and response characteristics. Results show that the weights of UECIIMP fluctuate slightly (within 10%) before and after sensitivity analysis, with good stability and reliability. UECIIMP in Hangzhou > Hefei ≈ Nanjing, spring ≈ autumn > summer ≫ winter. From 2009 to 2019, UECIIMP has improved in all 33 districts of the Hefei-Nanjing-Hangzhou Region. The significant improvement of UECIIMP in 2014–2019 is 4.3 times than that in 2009–2014. The correlation between UECIIMP and economic index indicates that economic development has a positive impact on the urban ecological environment. The significant degradation of UECIIMP in the urban expansion area demonstrates a negative impact on the local environment from urban expansion.

Instability of remote sensing ecological index and its optimisation for time frequency and scale

Evaluating regional ecological environmental quality (EEQ) using remote sensing is important for achieving sustainable development. The remote sensing ecological index (RSEI) is a typical EEQ evaluation model used to comprehensively reflect regional ecological quality. However, owing to its complete dependence on remote sensing image information, the RSEI also has inherent issues, including unstable time series and inconsistent resolutions of its four sub-indices. To address these problems, this study used the Google Earth Engine (GEE) to propose an optimisation notion for the harmonic analysis of time series (HANTS) coupled random forest (RF) model and determined the accuracy and image quality of the optimised RSEI (RSEIo) for the Yangtze River Basin (YRB). Based on the findings, HANTS could fill the gaps in the images, effectively reducing the noise and discrete degrees of the four indices and thereby improving the stability of the RSEI. The correlation coefficient (R) between the RSEI and RSEIo was 0.93, with a root mean square error (RMSE) of 0.13, indicating that RSEIo is reliable. The image quality assessment (based on contrast and information entropy) indicates that HANTS combined with the RF model can produce RSEIo images with higher definitions and richer information at a constant spatial resolution. Further, the pixel-by-pixel coefficient of variation evaluation indicates that the RSEIo image was highly stabilised, yielding higher numbers of effective RSEIo images without changing the temporal resolution. Compared with traditional RSEI calculations, the optimisation proposal herein could highlight ecological differences caused by topographic changes in the YRB, which would produce an RSEI closer to actual surface conditions. Further, this proposed method could be used to obtain more detailed ecological information at a constant spatiotemporal resolution, thereby meeting the needs of long-term ecological monitoring in large-scale regions.

Time-frequency optimization of RSEI: A case study of Yangtze River Basin

Remote Sensing Ecological Index (RSEI) is one of the most widely used ecological quality assessment indicators. Due to the noise caused by adverse atmoshperic conditions and other factors, the RSEI calculated from the original image usually has the phenomena of lack of information and unstable image quality. Therefore, based on Google Earth Engine (GEE) cloud platform, this study adopts three common data reconstruction algorithms firstly, namely: Savitory-Golay filter (SG), harmonic analysis of time series (HANTS), Whittaker Smoother (WS), which are used to reconstruct the original MODIS time series data in the Yangtze River Basin (YRB) from 2000 to 2020, in order to optimize the calculation process of RSEI. At the same time, three indicators (correlation coefficient (R), standard deviation (STD), root mean square error (RMSE)) are used for the accuracy evaluation. The results show that data reconstruction can fill gaps in RSEI, the reconstruction performance of WS and SG for four parameters is better than HANTS, and the four SG reconstructed sequences have the strongest correlation with the original sequences (R between 0.8 ∼ 1), while the WS reconstruction sequence has the lowest error value (both STD and RMSE are less than 1), both of them can correct the pixel value, which is conducive to maintaining the stability of RSEI in the temporal dimension; the RSEI produced by HANTS has the best accuracy, that is, R, STD, RMSE are respectively 0.898, 0.130, 0.104. As shown by the research, it is necessary to de-noise each parameter before synthesizing RSEI. This study can provide a theoretical basis for applying time-frequency algorithms to optimize the ecological monitoring performance of RSEI.

Integrating Remote Sensing and Meteorological Data to Predict Wheat Stripe Rust

Wheat stripe rust poses a serious threat to wheat production. An effective prediction method is important for food security. In this study, we developed a prediction model for wheat stripe rust based on vegetation indices and meteorological features. First, based on time-series Sentinel-2 remote sensing images and meteorological data, wheat phenology (jointing date) was estimated using the harmonic analysis of time-series combined with average cumulative temperature. Then, vegetation indices were extracted based on phenological information. Meteorological features were screened using correlation analysis combined with independent t-test analysis. Finally, a random forest (RF) was used to construct a prediction model for wheat stripe rust. The results showed that the RF model using the input combination (phenological information-based vegetation indices and meteorological features) produced a higher prediction accuracy and a kappa coefficient of 88.7% and 0.772, respectively. The prediction model using phenological information-based vegetation indices outperformed the prediction model using single-date image-based vegetation indices, and the overall accuracy improved from 62.9% to 78.4%. These results indicated that the method combining phenological information-based vegetation indices and meteorological features can be used for wheat stripe rust prediction. The results of the prediction model can provide guidance and suggestions for disease prevention in the study area.

Accurate Monitoring of Submerged Aquatic Vegetation in a Macrophytic Lake Using Time-Series Sentinel-2 Images

Submerged aquatic vegetation (SAV) is one of the most important biological groups in shallow lakes ecosystems, and it plays a vital role in stabilizing the structure and function of water ecosystems. The study area of this research is Baiyangdian, which is a typical macrophytic lake with complex land cover types. This research aims to solve the low accuracy problem of the remote sensing extraction of SAV, which is mainly caused by water level fluctuations, differences in life-history characteristics, and mixed-pixel phenomena. Here, we developed a phenology–pixel method to determine the spatial distribution of SAV and the start and end dates of its growing season by using all Sentinel-2 images collected over a year on the Google Earth Engine platform. The experimental results show the following: (1) The phenology–pixel algorithm can effectively identify the maximum spatial distribution and growth period of submerged aquatic vegetation in Baiyangdian Lake throughout the year. The unique normalized difference vegetation index (NDVI) peak characteristics of Potamogeton crispus from March to May were used to effectively distinguish it from the low Phragmites australis population. Textural features based on the modified normalized difference water index (MNDWI) index effectively removed the mixed-pixel phenomenon of macrophytic lakes (such as dikes and sparse reeds). (2) A complete five-day interval NDVI time-series dataset was obtained, which removes potential noise on the temporal scale and fills in noisy observations by the harmonic analysis of time series (HANTS) method. We determined the two phenological periods of typical SAV by analyzing the intrayear variation characteristics of NDVI and MNDWI. (3) Using field-survey data for accuracy verification, the overall accuracy of our method was determined to be 94.8%, and the user’s accuracy and producer’s accuracy were 93.3% and 87.3%, respectively. Determining the temporal and spatial distribution of different SAV populations provides important technical support for actively promoting the maintenance and reconstruction of lake and reservoir ecosystems.

A novel regression method for harmonic analysis of time series

Harmonic analysis of time series is an important technique to reveal seasonal land surface dynamics using remote sensing information. However, frequency selection in the harmonic analysis is often difficult because high-frequency components are useful for delineating seasonal dynamics but sensitive to noise and gaps in time series. On the other hand, it is challenging to obtain temporally continuous satellite data with high quality because of atmospheric contamination. We developed a novel regression method named Harmonic Adaptive Penalty Operator (HAPO) for harmonic analysis of unevenly distributed time series. We introduced a new penalty function to minimize unexpected fluctuations in the model, which can substantially reduce the overfitting issue of regression in time series with temporal gaps. Specifically, the new penalty function minimizes the length of the model curve and the value range difference between the model and time series observations. We compared HAPO with three widely used regression methods (OLS: Ordinary Least Squares; LASSO: Least Absolute Shrinkage and Selection Operator; and Ridge) with different scenarios using Landsat time series data across the United States. First, we evaluated methods using Landsat surface reflectance time series within a single year. HAPO showed small and consistent monthly Root Mean Square Deviation (RMSD) values, in which most of the time RMSD values of predicted reflectance were less than 0.04. More importantly, HAPO showed consistent and less bias given varying density and irregularity of time series. Second, we evaluated methods using multi-year time series and the result suggested that HAPO was a better predictor of relatively short time series (less than4 years) with steady small RMSD values. When a longer time series (≥4 years) was used, all four methods disclosed similar RMSD values, but HAPO outperformed other three methods when there were temporal gaps. Last, we preliminarily tested how regression methods affected change detection and classification accuracy. HAPO showed the highest change detection accuracy of all tests in terms of F1 score when using the change threshold of 0.9999. In classification, HAPO produced the highest accuracy for short time series segments (one- or two-year time series). In contrast, all methods reached similar accuracy for 5-year time series. These results suggest that for areas that have large seasonal observation gaps or for time series that have less than 4 years records, HAPO can provide more consistent and accurate analytical results than other regression methods for harmonic analysis of time series.

Comparison of Multi-Methods for Identifying Maize Phenology Using PhenoCams

Accurately identifying the phenology of summer maize is crucial for both cultivar breeding and fertilizer controlling in precision agriculture. In this study, daily RGB images covering the entire growth of summer maize were collected using phenocams at sites in Shangqiu (2018, 2019 and 2020) and Nanpi (2020) in China. Four phenological dates, including six leaves, booting, heading and maturity of summer maize, were pre-defined and extracted from the phenocam-based images. The spectral indices, textural indices and integrated spectral and textural indices were calculated using the improved adaptive feature-weighting method. The double logistic function, harmonic analysis of time series, Savitzky–Golay and spline interpolation were applied to filter these indices and pre-defined phenology was identified and compared with the ground observations. The results show that the DLF achieved the highest accuracy, with the coefficient of determination (R2) and the root-mean-square error (RMSE) being 0.86 and 9.32 days, respectively. The new index performed better than the single usage of spectral and textural indices, of which the R2 and RMSE were 0.92 and 9.38 days, respectively. The phenological extraction using the new index and double logistic function based on the PhenoCam data was effective and convenient, obtaining high accuracy. Therefore, it is recommended the adoption of the new index by integrating the spectral and textural indices for extracting maize phenology using PhenoCam data.

Bayesian Temporal Tensor Factorization-Based Interpolation for Time-Series Remote Sensing Data With Large-Area Missing Observations

Land surface temperature (LST) is widely used in the field of time-series remote sensing. However, due to the influence of cloud cover, the large area of LST data observation is missing, which seriously affects the later data analysis. In the past research, various effective interpolation methods have been developed, but they usually cannot effectively interpolate the image data with large observation missing. In this article, a new method for interpolating these missing data called Hilbert tensor rearrangement with Bayesian temporal tensor factorization (HTR-BTTF) is proposed. This method requires tensor rearrangement of remote sensing data, combined with BTTF method for interpolation. In order to evaluate the performance of our method, we select three real study areas with different climates, Wuhan, Harbin, and Kunming LST data during day and night, and add cloud covers with different sizes to the cold and warm season layers each year. BTTF, inverse distance weighted (IDW), harmonic analysis of time series (HANTS), and GapFill are used as comparison methods. Root-mean-square error (RMSE) is a comprehensive evaluation index of interpolation results. Experiments have shown that HTR-BTTF is an effective method for interpolating missing observations, which is better than other methods. In the simulation experiment of the largest cloud cover size, on average, the RMSE of the data filled using the HTR-BTTF method was 17.2&#x0025; lower than that of BTTF and 54.9&#x0025; lower than that of the GapFill method, and it shows good robustness and high accuracy.

Use of a BP Neural Network and Meteorological Data for Generating Spatiotemporally Continuous LAI Time Series

Spatiotemporally continuous long-term leaf area index (LAI) products are urgently needed to monitor environmental changes. The current filter- or curve-fitting-based time series reconstructive algorithms fail to reconstruct the LAI time series with many continuous missing values or missing values in key phenological periods, which are common issues in high-spatial-resolution LAI time series. In this article, a meteorological data-driven backpropagation neural network (MBPNN) was proposed to reconstruct discontinuous LAI profiles with a two-step process using vegetation phenological information. As the basis of the strong dependence of vegetation growth on meteorological conditions, a reasonable growth trajectory of reconstructed LAI can be guaranteed by the algorithm even though if many observed values are missing. Validations for reconstructed LAI were conducted both spatially and temporally based on reference maps and field-measured long-term observations. The results showed that the LAI predicted by the MBPNN had a similar accuracy (RMSE = 0.4076) as the Landsat LAI inversions (RMSE = 0.4083) and a similar reconstructed trajectory as the field-measured LAI series even though over 100 days of continuous data were missing (RMSE = 0.1620). A comparison with the Harmonic ANalysis of Time Series (HANTS) algorithm showed that the accuracy of MBPNN was more stable regardless of the size/position of the missing data, and the proposed method performed much better when the data were continuously missing for 50 days or more.

Optimal Estimate of Global Biome—Specific Parameter Settings to Reconstruct NDVI Time Series with the Harmonic ANalysis of Time Series (HANTS) Method

Terrestrial remote sensing data products retrieved from radiometric measurements in the optical and thermal infrared spectrum such as vegetation spectral indices can be heavily contaminated by atmospheric conditions, including cloud and aerosol layers. This contamination results in gaps or noisy observations. The harmonic analysis of time series (HANTS) has been widely used for time series reconstruction of remote sensing imagery in recent decades. To use HANTS model, a series of parameters, such as number of frequencies (NF), fitting error tolerance (FET), degree of over-determinedness (DoD), and regularization factor (Delta), need to be defined by users. These parameters provide flexibilities, but also make it difficult for non-expert users to determine appropriate settings for specific applications. This study systematically evaluated the reconstruction performance of the model under different parameter setting scenarios by simulating pixel-wise reference and noisy NDVI time series. The results of these numerical experiments were further used to identify optimal settings and improve global NDVI reconstruction performance. The results suggested optimal settings for different areas (local optimization). If a user opts to use unique settings for global reconstruction, the setting NF = 4, FET = 0.05, DoD = 5, and Delta = 0.5 can produce the best performance across all setting scenarios (global optimization). In addition, several internal improvements, such as dynamic weighting scheme, polynomial and inter-annual harmonic components, and ancillary attributes of input data can be used to further improve the performance of reconstruction. With these results, future non-expert users can easily determine appropriate settings of HANTS for specific applications in different regions.

Performance Evaluation of Long NDVI Timeseries from AVHRR, MODIS and Landsat Sensors over Landslide-Prone Locations in Qinghai-Tibetan Plateau

The existence of several NDVI products in Qinghai-Tibetan Plateau (QTP) makes it challenging to identify the ideal sensor for vegetation monitoring as an important factor for landslide detection studies. A pixel-based analysis of the NDVI time series was carried out to compare the performances of five NDVI products, including ETM+, OLI, MODIS Series, and AVHRR sensors in QTP. Harmonic analysis of time series and wavelet threshold denoising were used for reconstruction and denoising of the five NDVI datasets. Each sensor performance was assessed based on the behavioral similarity between the original and denoised NDVI time series, considering the preservation of the original shape and time series values by computing correlation coefficient (CC), mean absolute error (MAE), root mean square error (RMSE), and signal to noise ratio (SNR). Results indicated that the OLI slightly outperformed the other sensors in all performance metrics, especially in mosaic natural vegetation, grassland, and cropland, providing 0.973, 0.015, 0.022, and 27.220 in CC, MAE, RMSE, and SNR, respectively. AVHRR showed similar results to OLI, with the best results in the predominant type of land covers (needle-leaved, evergreen, closed to open). The MODIS series performs lower across all vegetation classes than the other sensors, which might be related to the higher number of artifacts observed in the original data. In addition to the satellite sensor comparison, the proposed analysis demonstrated the effectiveness and reliability of the implemented methodology for reconstructing and denoising different NDVI time series, indicating its suitability for long-term trend analysis of different natural land cover classes, vegetation monitoring, and change detection.